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Abs t rac t 


In the absence of an applied magnetic field a process which 

will permit the existence of a Fermi contact hyperfine interaction 
6 3 

at the Cu p nucleus which does not average to zero when summed 
over the randomly oriented spins of N conduction electrons is the 
scattering of s-wave electrons from the long-ranged exchange poten- 
tial of an open-shell impurity atom located at a distance R from 
the Cu nucleus and then the Fermi contact interaction at the Cu 
nucleus. This is a two-center analog of the virtual process which 
is responsible for most of the Fermi contact hyperfine structure 
of open-shell atoms (polarization of the closed s core). The first 
order energy due to the Fermi contact interaction of the conduction 
electrons at the Cu nucleus which have first exchange scattered 
from the 3d electrons of Mn is calculated as a function of R. 

It is found to change sign between R*l*.8299 (nearest neighbor 
distance), R»4.8299 and R=8, and R=8 and R«10 au, suggesting 
dependence on impurity concentration experiments and experiments 
in which compression or expansion can change the lattice distances. 
The exchange potential is presented as a function of electron 
velocity and R. For completeness the RKKY type spin density 
waves are presented as a function of the distance from the impurity 
for several values of R. Also for R=0 an estimate is made of the 
exchange polarization potential due to the adiabatic polarization 
of the impurity by the scattering electrons and found to be signi- 
ficant. 
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for the spin up electron and the absence of these potentials 

for the spin down electron. This is analogous to the virtual 

process which gives the major contribution to the Fermi contact 

hyperfine structure in open shell atoms due to the exchange of 

an s electron with an open shell electron, its scattering into 

an s excited state, then the contact interaction and the scattering 

5 6 

back into the original state. ’ This process is represented 
d i agrammat i ca 1 ly by, 



Therefore many conduction electrons with randomly oriented spins 

will experience a net nonzero Fermi interaction with a Cu nucleus 

in the vicinity of an open shell impurity atom, in the absence of 

an applied magnetic field. The strength of this interaction will 

♦ 

depend on the magnitude of the exchange potential and thus on the 
distance between the Cu and impurity atoms. The interaction energy 
can be of either sign and thus does not have to be in excess of 
the no-impurity shift. The connection with the HWJ6 and other 
theories now becomes more apparent. The above describes the micro- 
system of a single Cu-impurity pair. Since the other impurities 
of the system can be distributed among 2S + 1 degenerate states, 
the Fermi interactions over the entire system will still average 
to zero unless there is an applied field or short ranged magnetic 
order among the impurity atoms (known to exist for Mn^) . This 
paper will be concerned with a given pair for which the impurity 
has a particular spin projection (Mg*S is the case chosen) and 
thus with the zero-field shift localized in the region of this pair. 



(u?:er p a g a 



It is clear what happens when the field Is turned on. At zero 
temperature all of the Mn atoms occupy the lowest of the 2S+1 
Zeeman-split states, and the local shifts due to the Fermi 
Interaction are of the same sign and add. As the temperature 
increases the levels become equally populated, and the locai 
shifts average to zero over the entire system. This is in 
agreement with the observations of Sugawara. 


/ 
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Theo ry 


The first order energy due to the Fermi contact interaction 
at the Cu nucleus is given, by, 


£§< 4* ka (£.) |2M e ii N 8ir/36( < r)s,*i|^ k0 ( < r)> 


where represents the set of quantum numbers k&m (onlyil-m-0 
contribute to this integral by the operation of the delta 
f unct 1 on6 ( £) ) and a is an index over the spin states of N 
conduction electrons, each of which has a wave function 
(dropping the spin indices), 


4 k (r.) - (2/R o ) i kJa Jl (i) J! '(2Jl+l) i e l6 Jlu k>)l (r)P jl (cos0)/(4Tr) i 2 


where the radial waves are normalized by tie condition that they 
vanish at the boundary of a sphere of radius R^; hence. 


i - V ,/dk i. 
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Is the phase shift, defined as usual by the condition, 


lim u^ ^(r) " sln(kr - tt/ 2A + 6^)/kr 
r->“ * 


and a^ are defined such that a^a^ is the Fermi distribution 
function which at zero temperature (assumed throughout) if 1 
for electrons whose energies are below and 0 for electrons 
whose energies are above, where e is the Fermi energy, taken 


to be 7ev for the free electron model of Cu 
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s and t are the 


electron and nuclear spin angular momenta respectively, 

y.. are eti/2m c and 2.22617 efi/2m c respectively, and 
n e p 


and 


6(0 - g6(r)/r 2 (2A+1)P A (cos0)/^Tr 


The particle wave function is given in the first Born iteration 

by , 


* k <r> - .'Vi - /dr's(r.r')V(i ).'Vi 
6(r,r ) - ^ 

I r-r 1 | 


6 
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The potential is defined, 


W,t(N*."0 A 

V(r) - 2(-2 1 /r Mn - Z 2 /r Cu ♦ £/d r' ■ . ( r' ■ ) ( 1 ±6 (m % , m g ' ) ) P , ^.(r’ ) 


r - r' 


* ^ (n» s »n» s ' )P { j'i'j (£.' ) (ej - k 2 /2)>) 


where m and m' are the spin projections of the scattering and 
atomic electrons respectively, Z ^ , N ^ and Z^, are the charges 
and numbers of electrons of Mn and Cu respectively, £j are the 
orbital energies of the bound $tatesij/j (_r) on Cu and Mn, and Pjj 
permutes the coordinates of the bound and scattering electrons, 
giving exchange potentials for (+) singlet and (-) triplet 
scattering respectively. The key point of the potential is the 
fact that when the summation over spin states is performed in 1 
there is cancellation of the products of the Fermi interaction 
and the coulomb (direct) interaction and cancellation of the 
products of the Fermi interaction and the exchange interaction 
for all interactions except those products which involve exchange 
with the Mn d electrons. 

Rigorously we should consider the scattering as modified 
by the residual l/ r £ u (coulomb tail) potential left when the ks 
electron is removed and the band structure (resulting from use 
of Bloch waves) which the less completely screened potential can 
support at close distances. The first consideration is not relevant 
in the first Born method of solution used here (iteration with un~ 
distorted plane waves), and the neglect of band effects is at least 
qualitatively correct for the long-ranged exchange process considered 
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here. Further there Is an additional potential, V j ± ^pol exch' 
which results from the long-ranged polarization of the target by 
the scattering electron and is known to affect the phase shift by 
about \0% In the low energy elastic scattering from hydrogen atoms. 
We have estimated the contribution of V j exc h for t * ie fictitious 
case (R-0) when the Cu and Nn atoms are united and found the 
exchange potential changed by about 10$. A careful study of this 
potential at nonzero R is therefore Indicated (see the appendix) 

We specialize eq. 6 to r>0 (the only contribution to the 
delta function interaction in 1), leaving only the leading term 
in the expansion of G, 


lim^.(r) -(1 - 1/lnr/dr'e Jkr /r 'VU') e *-*- ) x (2/R.) *ka ft e 1 5 0/ (4 tt) * 
r-f 0 - ” “ o u 
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where V. is obtained by keeping cnly the s-d exchange terms In 
7, since all others cancel when 1 is evaluated. We then define, 


J(k,R)/k - -l/Wd^'cosfkr'J/r'VU'JjgUr') 
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where we have kept the leading term in the plane wave expansion 

of e . This is equivalent to spherically averaging the potential 

so that the £,-0 partial wave, which is the only wave wh i ch 'cOn t r i • 
butes to 1, is uncoupled to other partial waves. 

If the spin up states are given by 8, then the spin down states 
are given by , 


lim (rj - 1x(2/R n )*ka n c ,5 0/(4Tr)* 
r*0 — oo 


10 
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which is a consequence of the choice Mg-S for the ground 
state of Mn. Thus to first order in J eq . 1 becomes. 


AE 


* F 

4/TTjJak k J(k,R) 


1 1 


when the exchange potentials take the upper sign (singlet 
scattering). The total shift is given by, 


AE - 1/I»AE S + 3AAE* 


12 


where AE C results when the potentials take the upper sign 
(triplet scattering). 
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. RKKV Sp t n Dens 1 ty 


The spin density function is defined here as p+~pt, 
where , 


pt + - |\|> k ++| 2 


I I 

■ e - e'* <r / 1 »irr / d£ e • — (V.(£, ) +2 (8Tr/3P e M N 26 (_r 
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for large distances from the Cu nucleus. For close distances we 
can write. 


* k (D 


ik r 
e 


I/inr/dr 


ikr /r (V ( r ) +2 


(8ir/3u e u N 26 (r_ 


) «n $ m | ) 


*k > 
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Now substitute 14 into the right hand side of 13 and keep only 
cross products of the Fermi and interactions. When the modulus 
is taken and the difference of spin down and spin up densities 
summed over k, the result to first order in J is, 


2 

o(r) » Vir^ak s l n (2kr) /r J (k , R) x4/3U e y N m s m| xl /Air 


r ■ (r Hn + 1(2 ' Rr Hn RcosR M^ 
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we average this expression over the angles of Mn in the usual 
way and define the perturbation to the electron density resulting 
from the two-center scattering process as a function of the 
distance from the impurity atom as, 


p ( r ) ■ 1/4ir/d^d0s in8o(/) 
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This function is plotted in Figs. 1, 2 and 3 for several 
values of R. 

it is important to note that even though eq . 13 results 
from the large r expansion of the Green's function (6), it 

I 

is the only surviving term in the limit r »0 when the Fermi 
interaction in 13 is evaluated, and likewise for the small 
r expansion of the Green'si function to give H. Thus the 
result ( 15 )! is not asymtotic for large r. 
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Numeri cal Procedures 


The exchange Integrals defined by j(k,R) are the two-electron, 
two-center type, and the integrations are performed numerically 
in prolate spheroidal coordinates defined by. 


« * <r «n +r C )/R 
n ' (r Mn' r Cu )/R 


where R is along the polar axis of the reference frame. The 
inverse distance } I are 9' ven spheroidal coordinates 
by the Neumann expansion^, 

2/REE(2j.1)((j-H>!/<.M«l>! )M"I ( ?< ) ft l«l ( S> ) p j«l ( n ( ) pi”) („ 2 ) 

x (-l) m 18 


where P and Q are the associated Legendre functions of the first 
and second kind respectively. The results were checked by taking 
R*.01 in the spheroidal coordinates and checking against the R=0 
integrals in spherical polar coordinates. Note that in the R«0 
limit only the j»2 term of the Neumann series is nonvanishing in 
the integral defined by 9 owing to the angular integrations over 
and s and d pair. For RfO, however, the expansion must be carried 
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to convergence , and the contribution for each term through 
j-8 is tabulated in Table 2 . The overlap terms in eq. 7 
were estimated to contribute about 1$ at R-1 to the energy 
shift and were neglected (note that only the m«0 d electron 
has an overlap with the s waves). 

Inspection of Table 2 shows that for R greater than about 
4 the convergence is very slow. The numbers for the higher R's 
are not fully converged, but are adequate to establish that the 
oscillatory variation of the shift with R is not an artifact of 
a nonconverged series. Also it is felt that this convergence is 
adequate given the set of physical approximations used in this 
calculation. 

All integrals were evaluated using the analytic Hartree- 
Fock d orbitals of Mn calculated by Clement!.^ 


Results and Discussion 


Values of J(k,R) and AE S are tabulated in Tables 1 and 2. 

Two interesting features of these results are the very long 
range of the exchange potential and the oscillatory dependence 
of the energy shift on R. The long range is familiar in the 
evaluation of atomic hyperfine structure by perturbation 
theory^'^ in which better than 90$ of the contribution to the 
core polarization exchange diagram (see the Introduction) comes 
from the continuum excited states of the atom. This is due, of 
course, to the unboundedness of the continuum functions; thus the 
boundedness of the integrals is governed by how tightly the d 
electrons are bound. The oscillatory behavior of AE S can be 
inferred from the oscillatory dependence of J(k,R) on k for 
larger R by inspection of Table 1, and is to be expected when 
the interaction is small at large R, leaving a kind of distorted 
overlap of d orbitals with continuum orbitals. This means that 
there will be a complicated dependence of the shift on impurity 
concentration and suggests dependence on concentration experiments 
and/or experiments in which compression or expansion of the metal 
can change the lattice distances. Note that the interaction is 
almost at a node for the nearest neighbor distance R«A.8299 and 
is larger and of opposite sign at the next nearest neighbor 
distance R-6.8305. 

In the interpretation of the results for larger R it Is 
important to appreciate the smallness of the interaction and 
thus the tentativeness of the numbers of this calculation, from 
the point of view of the absolute numerical accuracy but especially 
from the point of view that better approximations to the physics 
need to be considered. Some of these include the exact numerical 
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solution of the relevant equations to obtain scattering wave 
functions more accurate than the first Born iteration used 
here, the consideration of band effects, the adiabatic polar- 
ization of the target (see the appendix), electron correlation, 
and the deviation of the d orbitals from those of the free 
atom. Also in Cu-Mn it is very important to consider the 
short ranged magnetic order (clustering) of the solute atoms. 
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.002826 -.003^57 .003119 -.002141 



Table 2. -AE /\x for selected R(a.u«). j Is the index In the expansion of 
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Figure Captions 


Fig. 1. p(r)/a where a- 8S 8 function of r, the 
distance from Mn, averaged over the angles of Mn, for several 
values of R. Calculated using d vh 1 ues ..conve rged t th rough ’ j *2 . 

Fig. .2. p(r)/a where a- ^x.j/5 as a function of r, the distance 

from Mn, averaged over the angles of Mn, for several values of 
R. Calculated using J values converged through j**2. 


Fig. 3* p(r)/a where as a function of r, the distance 

from Mn, averaged over the angles of Mn, for several calues of 
R. Calculated using J values converged through j«2. 
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Appendix. Adiabatic exchange polarization potential. 

The contribution of the adiabatic exchange polarization 
potential is estimated as follows. When the Mn target is 
polarized by a static electron at r^, the solution to the 
first order equation from Ray lei gh-Schrodi nger per turba t ion 
theory ; (for a perturbation on the m-0 orbital in particular, 
where nondegenerate perturbation theory can be used because 
the d orbitals are not split by the dipole potential of 

Q 

polarized orb i t.a 1 theory 7 ) is, 

• • -(2) 7/2 /(72O) i x(5/^nr) i u 1 (r)/rP 1 (?*? 2 )/r^P 2 (cos0)e(r<r 2 ) 1 


wheree is a step function 1 for r<r 2 and 0 for. r>r 2 where the 
radial part of the 3 d orbital (unnormal i zed) i s taken to 
be r exp(-r) (based on a binding energy of about .6 a.u. as 
calculated by Clementi). We do not take dementi's HF function 
for the unperturbed orbital because then we c’ould not find 
an analytic solution to the first. order equation. The first 
order radial solution is, 


u t (r) - e’ r (r 2 +r 3 /2+rV5*r 5 /1.5+11i l f;|lHx15) + 1 1.r 7 /(4xl4x15) 

+22 r 8 /(9x4xl4x1 5) +154r 9 /( 35x9x4x1 4x1 5>+308r 10 /(l 1x35x9x4x14x15) 


+ 0(10 -i *) r 11 ) 


2 
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The exchange polarization potential is generated by substituing 
1 for one tti ^ in eq . 7 of the text. J(k,R«0) is about 102 of 
j(k,R*0) calculated with unpolartzed orbitals. It is clear 
that this effect should be looked at more quantitatively. 






